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Motivated by the results of two-dimensional conformal field theory (CFT) we investigate the finite-size scaling 
of the mass spectrum of an Ising model on three-dimensional lattices with a spherical cross section. Using a 
cluster-update Monte Carlo technique we find a linear relation between the masses and the corresponding scaling 
dimensions, in complete analogy to the situation in two dimensions. Amplitude ratios as well as the amplitudes 
themselves appear to be universal in this case. 



1. INTRODUCTION 

Divergent correlation lengths of a statistical 
system approaching a critical point establish the 
symmetry of scale invariance on which finite-size 
scaling theory is based. In a continuum critical 
theory, however, there are additional symmetries 
present, namely translational, rotational and in- 
version invariance. The implications of the group 
of these symmetries, as explored in conformal 
field theory [|l]j|], go far beyond simple scaling 
theory, especially for two-dimensional systems, 
where the whole content of scaling operators as 
well as for example higher correlation functions 
can be generically extracted from this symmetry. 

As for a prominent result, conformal invariance 
suffices to determine the critical two-point func- 
tion in the plane. Using a conformal (logarithmic) 
transformation that wraps the plane around an 
infinite-length cylinder one thus arrives at an ex- 
pression for the correlation function of a primary 
(conformally covariant) operator on the cylinder 
|H , a geometry one generically considers in trans- 
fer matrix calculations. The fact that this trans- 
formation does only affect the radial but not the 
angular part of the coordinates led Cardy H to a 
generalization for the higher-dimensional geome- 
tries S 1 x IR. Thus, exactly as in two dimen- 
sions, the large distance exponential decay of cor- 
relations is determined by a correlation length 
that is given by the inverse of the scaling dimen- 



sion x of the operator under consideration, 
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where R is the radius of S . It has to be noted, 
however, that the meaning of a primary operator 
is a priori not well defined in three dimensionsn so 
that this result has to be considered a conjecture 
and a numerical analysis is more than an exercise 
in this case. 

A first attempt to establish such evidence using 
the Hamiltonian formulation of the Ising model 
and Platonic solids as discretization of the sphere 
S 2 was inconclusive due to the restricted num- 
ber and size of regular polyhedra H . Using ap- 
proximate quadrangulations we are able to check 
Cardy's conjecture for systems sufficiently large 
to carry out a proper finite-size scaling (FSS) 
analysis. 

2. LATTICE DISCRETIZATION 

There are several straightforward choices for 
model lattices with spherical topology J7|; the 
more natural are variants of a rectangular mesh 
on a cube. Different refinements can be applied to 
reduce the concentration of the curvature around 
the corners of the cube such as filling in triangles 
instead of the cube corners or projecting the cube 
on the sphere by varying the link lengths, i.e., ap- 
plying an appropriate site-dependent weight func- 
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therein. 
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tion to the action. For those cube-like lattices, 
however, differences in the scaling behavior of 
bulk quantities are found to be quite small |j. 
As for our observables, there is some evidence to 
believe that ratios of correlation lengths of pri- 
mary operators are universal []S|-[10|], so that we 
can expect good agreement regardless of the lat- 
tice used if Cardy's conjecture holds. For the 
amplitudes themselves one might find effects of 
the non-uniformity of the lattice, but would not 
expect them to become visible until a very high 
level of precision is reached. We thus here use the 
simple cubic approximation to the sphere by six 
Lx L square lattices with suitably chosen bound- 
ary conditions; effects of the discretization will be 
considered elsewhere UM. 

Having chosen a lattice approximation to the 
sphere one has to fix a definition of the radius R 
of the sphere the cube with edge length L should 
correspond to. As is easily checked the three pos- 
sibilities of assigning a unit area to each lattice 
site, each pair of bonds, or each lattice square, 
respectively, lead to total lattice areas of 

, 6L(L-2)+8 "sites", 

1 6L(L-2)+6 "bonds", "squares", [) 



and generate via the relation R = y^A/An two 
different sorts of pseudo radii, which only differ by 
a constant shift, thus leading to slightly different 
approaches to the leading FSS amplitude. 

3. SIMULATIONAL DETAILS 

We consider a classical, ferromagnetic, nearest- 
neighbor Ising model with Hamiltonian 



H 



-J 2_^ SiSj, 
(i,3) 
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on lattices compound of the above described cu- 
bic discretizations of the sphere times a linear di- 
rection of length L z , modelling the K-direction. 
In order to minimize the effect of a finite L z it 
was (self-consistently) chosen so that L z /£ « 15 
and periodic boundary conditions in z-direction 
were applied. Simulations were done using the 
Wolff single-cluster update algorithm at an in- 
verse temperature (3 = 0.2216544(3) fTJ], check- 
ing for the influence of a temperature shift off the 



critical point by reweighting. To enable a proper 
FSS analysis system sizes ranging from L = 4 to 
L = 12 were used, corresponding to radii from 
about 2 to 8 and overall lattice volumes of 6 000 
to 300 000 spins. 

Since the densities of energy and magnetiza- 
tion are (the only) primary operators in the two- 
dimensional Ising model and the lowest lying 
states of the two sectors, we here consider the am- 
plitudes of the corresponding (exponential) cor- 
relation lengths A a / e = ^a/t/R and their ra- 
tio. Simulational measurements were done for the 
(connected) correlation functions G c (z) of these 
operators using a zero mode projection. Estima- 
tors for the correlation lengths that eliminate ad- 
ditive and multiplicative constants in the expo- 
nential decay are then given by 

G c (z)-G c (z-Ay~~' 
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G c {z + A) - G c {z) 

with A = 1, 2, 3, . . ., from which overall estimates 
£,(R) are formed by averaging. A jackknife tech- 
nique was used to reduce estimator bias and for 
the error analysis; for details see Jl0| , p^[ . 

4. NUMERICAL RESULTS 

Using the above described procedures one ends 
up with scaling plots for the correlation lengths 
as shown in Fig.|lj(a) for the magnetization den- 
sity that show an almost perfect linear behav- 
ior and invisible differences between the different 
radii definitions according to Eq. (g). Plotting the 
amplitudes £/i? (Fig. 0(b)), however, reveals the 
presence of corrections to scaling that enforce the 
use of nonlinear fits of the form 



£(i?) =AR + BR a , 



(5) 



with a free parameter a that introduces a sys- 
tematic error to the final amplitudes A a / C due to 
higher order corrections. A combination of the 
fits for the two sorts of radii gives the amplitudes 



A a = 1.996(20), 
A e = 0.710(38), 
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which agree well with the conjectured amplitudes 
A conj = 1 j Xg = 1.9324(19), 
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0.70711(35), 
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Figure 1. (a) FSS plot for the spin correlation 
length £,a(R)- (b) Scaling of the amplitudes t; a /R- 
The horizontal line indicates the conjectured am- 
plitude according to Eq. (^). 



taking into account the systematic error due to 
the different radii definitions that is of the same 
order of magnitude as the statistical error. The 
amplitude ratios compare as: 



A a /A e 
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2.81(15), 
2.7326(16). 



(8) 



5. CONCLUSIONS 



The presumably universal ratio of the ampli- 
tudes of the spin and energy correlation lengths of 
the Ising model on a cubic model of S 2 x IR agrees 
with the inverse ratio of the scaling dimensions 
as conjectured by Cardy for the continuum case. 
Moreover, corrections to scaling due to the non- 



uniform distribution of curvature over the lattice 
that might influence the amplitudes themselves 
seem not to be very important at the given level 
of accuracy; only the slightly shifted result for 
the spin amplitude A a might indicate the onset 
of such effects. In connection with the results 
for systems with toroidal cross section ]H],|l3| this 
seems to indicate a deeper analogy between the 
2D and 3D situations. In this latter case, how- 
ever, the system is not conformally related to a 
flat space which renders impossible an analytical 
treatment of Cardy's kind. In the spherical case, 
on the other hand, conformal flatness is fulfilled, 
but the algebraic meaning of the primarity of an 
operator remains unclear. Thus our result shows 
that spin and energy densities of the 3D Ising 
model are primary operators, taking Eq. (Q) as a 
definition of primarity in three dimensions. 
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